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I. INTRODUCTION 



A fundamental goal of perturbative QCD is to provide a microscopic description of 
hadronic short-distance phenomena that yields reliable predictions to be compared with 
experimental data of increasing precision. While singularities on the timelike axis in the 
complex Q 2 plane of hadronic observables are related to physical particles (or resonances), 
the appearance of singularities on the spacelike axis are unphysical and may violate causal- 
ity. On the other hand, the expansion of hadronic quantities at large momentum transfer Q 2 
can be safely calculated in terms of a power-series expansion in the running strong coupling 
a s {Q 2 ) by virtue of asymptotic freedom. But the one-loop running coupling contains at 
Q 2 = Aq CD (Aqcd = A in the following) a ghost singularity — the Landau pole — that spoils 
its analyticity structure. To restore analyticity and ensure causality in the whole Q 2 plane, 
this pole has to be removed. With most available experimental data on several exclusive 
processes being at rather low Q 2 values, the Landau-singularity problem is not only of aca- 
demic interest, but affects significantly perturbative predictions in the low-to- medium Q 2 
domain. The reason is that — lacking all-order perturbative expressions — one has to resort 
to a renormalization-scheme choice that makes the uncalculated higher-order corrections 
negligible and adopt a renormalization scale that reflects the typical parton virtualities in 
the considered process. The latter procedure, however, may result into a scale in the region 
of only a few A, where the application of perturbation theory for the conventional running 
coupling without infrared (IR) protection against the Landau pole becomes inapplicable — a 
prominent example being the Brodsky-Lepage-Mackenzie scale-setting procedure [1]. Dif- 
ferent strategies have been suggested over the years as how to minimize the dependence 
on the renormalization scheme and scale setting — unavoidable in any perturbative calcula- 
tion beyond the leading order — and obtain reliable and stable results in the low-momentum 
regime (see, for example, Ref. [2] for a recent extensive discussion of these issues in terms of 
the electromagnetic pion form factor and references cited therein). 

In a series of papers during the last few years Shirkov and Solovtsov (SS) [3-7] have 
developed an approach which enables the removal of the Landau singularity without in- 
troducing extraneous IR regulators, like an effective gluon mass [8-12]. The analyticity of 
the coupling in the spacelike region is achieved by a nonperturbative, power-behaved term 
that contains no other scale than A and leaves the ultraviolet (UV) behavior of the running 
coupling unchanged. At zero-momentum transfer the Shirkov-Solovtsov coupling assumes 
a universal value that depends solely on renormalizat ion-group constants. Using dispersion 
relations, this scheme was both generalized (in approximate form) to higher-loop orders and 
also extended to the timelike regime [6, 7, 13-23], encompassing previous incomplete at- 
tempts [24, 25] in this direction, and amounting to the theoretical framework of Analytic 
Perturbation Theory (APT). There have been a number of parallel developments by various 
authors during the past several years to avoid the Landau pole using different "analytization" 
techniques, prime examples being Refs. [17-19, 24-35]. 

Two other major challenges, connected with — first — the implementation of the Shirkov- 
Solovtsov "analytization" to three-point functions beyond the leading order of perturbation 
theory and — second — the extension to non-integer powers of the coupling, remained open, or 
at least partially open. Indeed, in the first case, extensive analyses [2, 36, 37] have shown that 
the "analytization" principle has to be generalized to accommodate a second scale, serving 
as a factorization scale, or in order to include evolution effects comprising typical logarithms 
to some fractional power. Technically speaking, this means to extend the assertion of an- 
alyticity from the level of the coupling (and its powers) to the level of the whole reaction 
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amplitude. This requirement was formalized by Karanikas and Stefanis (KS) in [38, 39] in 
an attempt to calculate power corrections to the pion form factor and the Drell-Yan pro- 
cess. The systematic development of a perturbative expansion in terms of fractional powers 
of the coupling — the second major challenge — is the goal of the present investigation, the 
main focus being placed on the methodology towards improving perturbative higher-order 
calculations in QCD. This goal has been accomplished and will be described in this paper. 
A specific application of the KS "analytization" principle to the pion's electromagnetic form 
factor at NLO accuracy is given in fully worked out detail in [40]. Other applications will 
follow in future publications in conjunction with the inclusion of heavy-flavor thresholds and 
the extension to the timelike regime. 

The paper is organized as follows. In Sec. II we first review the key features of the 
original Analytic Perturbation Theory of Shirkov and Solovtsov, highlighting those proper- 
ties pertaining to the generalization of the approach to fractional powers of the coupling. 
The actual extension of the approach to fractional — in fact, real — powers of the coupling is 
performed in Sec. III. This section describes in three subsections the new "analytization" 
technique, based on the Laplace transform, the verification of the analytic properties of the 
obtained results, and the way to include products of powers of the coupling with powers of 
logarithms. Moreover, we provide here approximate expressions for two-loop quantities in 
terms of one-loop analytic-coupling images and their index derivatives that can be extremely 
useful in practical calculations. Section IV is devoted to the validation of the developed the- 
oretical framework of the Fractional Analytic Perturbation Theory (FAPT) and includes a 
table where we collect the algorithmic rules to connect the new analytic framework to the 
standard QCD perturbative power-series expansion. Our conclusions are drawn in Sec. V, 
while important technical details are presented in three appendices. 

II. ORIGINAL ANALYTIC PERTURBATION THEORY 

In the analytic perturbation-theory approach of Shirkov and Solovtsov, the power-series 
expansion in the running coupling is given up in favor of a non-power series (functional) 
expansion. This can be written generically in terms of numerical coefficients d m in the 
following way [4, 41] 



where the "normalized" coupling a = boa s / (4ir) (bo is the first coefficient of the QCD f3- 
function — see Appendix B) has been introduced instead of a s in order to simplify intermedi- 
ate calculations and because then the analytic coupling A\ is bounded from above by unity 
[4]. In the above expression, the superscript m on appears on the left-hand side (LHS) 

as a power, whereas on the right-hand side (RHS) the subscript m on Am denotes the index 
of the functional expansion; 1 (/) denotes the loop order. For the sake of simplicity, we will 
avoid to indicate the loop-order index explicitly because we mostly work in the one-loop 
approximation; deviations, if needed, will be labelled by appropriate superscripts or sub- 
scripts in parentheses, like in Eq. (2.1). The conversion to analytic images of the coupling 
is achieved in terms of the functions 




(2.1) 



Hi 



m 



an 



(2.2) 



In the following, a calligraphic notation is used to denote analytic images. 
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according to the general prescription 

^ )L = ip=M^U. (2 . 3) 

For the one-loop running coupling 



a ~ in(gVA 2 ) ' (2 ' 4) 



we have 



A(Q 2 ) - [.««)]„ - 1^X5) - qvX^T ^ 

and ^4^(0) = 1. Employing the variable L = ln(Q 2 /A 2 ), which naturally appears in pertur- 
bative QCD (pQCD) calculations, we can recast a and A in terms of L to obtain 

a 1 (L) = i A l {L) = j-^— l . (2.6) 

In this context, amplitudes (depending on a single scale Q 2 ) perturbatively expanded in 
terms of the powers of the running coupling map on a non-power series expansion [4, 41]: 

F(a) = J2 fna n (L) =► T{L) = £ f n A n (L) , (2.7) 

n n 

where /„ are numbers in minimal subtraction renormalization schemes. By construction, the 
set {A n } constitutes a linear space, which, however, is not equipped with the multiplication 
operation of its elements. Therefore, the product A n ■ A m has no rigorous meaning here. 
The standard algebra is recovered only for the main asymptotic contribution (cf. Eq. (2.6)) 
at L — > oo, when {A n } — > {a n } [3-5]. 

Let us now turn to the properties of this map and of the space {v4.„}. There are several 
points to note about them. 

1. The map should have the property of isomorphism, i. e., it should conserve the linear 
structure of the original space: 

a =>. A = 1 . (2.8) 

2. Renormalization-group summation leads to contributions like 

f(a) = a u , where v is real, (2.9) 

necessitating the introduction of the analytic images of f(a): [f(a)] an = [a"] an . These 
are exactly those terms needed to supply the original linear space {A n } with the 
completeness property as regards the differential operator with respect to the real 
index v. 

3. Motivated by the typical logarithmic contributions, appearing in loop calculations in 
standard pQCD, we consider the "analytization" of terms of the sort 

/«*)={ S^'Uu, < 2 - 10 > 

giving rise to the corresponding analytic images 

C v {a) = K ln(a)] a „, C v , m (L) ee KL m ] an . (2.11) 



4 



Expression (2.9) is universal and allows one to apply any one-loop renormalization-group 
results to APT. In fact, the corresponding renormalization factor Z, associated with the 
renormalizable quantity B, B{Q 2 ) = B (/i 2 ) Z \Q 2 ) / Z (/i 2 ) , reduces in the one-loop approxi- 
mation to 

Z ~ - 1(01, \ ' 

v = v Q = 70/ (26o) 

where 70 is the coefficient of the one-loop anomalous dimension. Therefore, we have 

[B(Q 2 )] an -[a^L)] an = A u (L). (2.12) 

The next two functions C u (a) and C u ^ m (L) appear in NLO of pQCD and also in light-cone 
sum rules [42, 43] and reflect the specific features of these calculations. An example of the 
first kind in connection with a NLO calculation of the electromagnetic pion form factor is 
treated in [40], while the investigation of such terms in the context of light-cone sum rules 
will be considered in a future publication. To be more specific, we will consider below terms 
of the form (a(2)) L and (a(2)) L 2 — rather than deal with the series containing the constant 
coefficients /„, given in Eq. (2.7). 

One possible way in generalizing the presented original APT formalism to non-integer 
(fractional) values of the index v is to construct the spectral density 

p u (a) = -Im [a u (-a)] (2.13) 



7T 



for i/£l. Indeed, substituting 

a( - a) = Tirrl ,V ' ^ = ln((T/A2) ' (2 - M) 

for the one- loop running coupling into Eq. (2.13), we can obtain by a straightforward calcu- 
lation a closed-form expression for the spectral density in the form 

. . 1 sm(u (p) ( L(a) \ /n . 

pM = ~ r 2 , T 2,w/2 ' v = arccos Trrhh 2 ■ ( 2 - 15 ) 

[tt 2 + L 2 (a)} 1 \\/L 2 {o) + 7r 2 y 

A result similar to that has been derived in the context of Electrodynamics in the early 
article of Ref. [44]. It was later re-invented in QCD by Oehme [45] and used by Magradze 
in [20]. To get now the desired analytic coupling for some fractional index, one has to insert 
this expression back into Eq. (2.3) and perform the integral numerically, loosing, alas, this 
way the possibility to reveal the mathematical properties of this function. Let us emphasize 
at this point that the extension of this procedure to the two-loop order for the first integer 
values of v has been done in Refs. [4, 5, 20], while the inclusion of still higher-loops [46] 
seems feasible. 2 However, this approach, based on the spectral density (2.13), is restricted 
to the specific structure of the Shirkov-Solovtsov APT. 



2 Indeed a partial result for a few fractional v values has already been obtained — Shirkov, private commu- 
nication. 
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III. FRACTIONAL ANALYTIC PERTURBATION THEORY 



A. A new generalization technique to include fractional indices 

In this subsection, we formulate and outline another procedure to continue the integer 
index of the analytic coupling to fractional values. First, to generate higher indices at the 
one-loop level of the analytic images within APT, or, equivalently, higher powers of the 
standard running coupling within conventional pQCD, we follow [4] and write 

ML)\_ 1 ( dY' 1 (Ai{L)\ 



a n (L) ) (n- 1)! \ dL J V al ( L ) 

Second, to facilitate the transition to fractional index values, it is instrumental to employ the 
Laplace representation of both types of couplings — the analytic, Ai(L), and the conventional 
one, a(L), (both at the one-loop order) — and define at L > 

^ff) = fW-^V (3.2) 
aHL)) y„ WO/ 

The advantage of this representation is that it transforms the result of a differential operator 
into an algebraic expression containing monomials. Then, applying Eq. (3.1) to (3.2), we 
get 



e 



t n-l 

(n- 1)! 



dt, (3.3) 



so that we establish the correspondence 



in—l 



A n (L) < > An(t) — ~, — — ttt • Ai(t) , (3.4) 

[n l). 

whereas for the case of the conventional pQCD coupling, one has the evident Laplace con- 
jugates a n 

a l {L) = - < — >■ fii = 1, (3.5) 

Lj 

j.n— 1 

a n {L) < — > a n = - — ■ ai . (3.6) 

Equation (3.4) enables us to generalize A n (L) to any real index v. To do so, let us 
introduce the following definition for the Laplace conjugate A u (t): 

Mt)= f ^yMt). (3.7) 

At this stage of the continuation in the index is, we have based our considerations solely on 
the first relation in Eq. (3.3). Therefore, the Laplace conjugate (3.7) remains valid for any 
non-power perturbative expansion satisfying this relation, reiterating that this holds true at 
the one-loop level. To complete the generalization process, we should obtain an expression 
for Ai(t), based on Eq. (2.6). This gives the result 

ML) = -- -j—j <— > Mt) = 1 - E S(t - m) , (3.8) 

m=l 
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which can be verified by a straightforward calculation. 

Let us pause for a moment to make some useful remarks concerning the behavior of 
the two parts of Eq. (3.8). One should note the strong difference in the behavior of these 
functions with respect to the logarithmic term of standard perturbation theory, on the one 
hand, 

and the pole remover appearing in APT, on the other, 

1 oo 
m=l 

Thus, one can define A V [V) according to Eq. (3.3), and, then, using Eqs. (3.7) and (3.8), 
arrive at 

/OO iu— 1 _ 1 1 00 

e~ Lt W) ■ Mt) dt = - - ^ • £ e- L - . (3.9) 

The series on the RHS of the latter equality coincides with the definition of the Lerch 
transcendental function [47] i/, i) at v' — v — 1 < for i — 1, i.e., 

E ^7 = ^(^ 1 -^ 1 )- ( 3 - 10 ) 

m=l 

The analytic continuation of s, 1) in the variables s, adopting the notation of Bate- 
mann and Erdelyi [47], determines $ as an analytic function of the variable z in the plane 
with a cut along (1, oo) for any fixed s (see for more details in Appendix A. 3 ). Finally, A v in 
Eq. (3.9) can be rewritten in the form of an analytic function with respect to both variables 
v and L; viz., 

A(£) = ^-^$(e- L ,l-M). (3.11) 
We state here and prove in Appendix A that A v is an entire function in v. 



B. Analytic properties 

To assess the analytic properties of Eq. (3.11), it is useful to recast the Lerch function 
vi a (see [47], Eq. (1.10.14) and also [48], Chapt. 8) 

z${z,v,\) = F{z,v) (3.12) 

entailing 



3 The transcendental Lerch function s, 1) is included in the widespread programs "Mathematica 5" and 
"Maple 7" . 
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where the first term in Eq. (3.13) corresponds to the standard PT, while the second one 
expresses the pole remover. Note that for a positive integer index, v = to > 2, one has the 
relation [47] 

F(z, 1 - to) = (-l) m F 1 - toJ , (3.14) 

so that substituting Eq. (3.14) in (3.13), one arrives at 

A m (L) = (-l) m A m (-L) (3.15) 

that confirms the specific symmetry relations worked out in [4]. From relation (3.15) and 
Eq. (A6) one obtains the explicit asymptotic expression for A m (L) at L — ► — oo 

AJL - -oo) = {-\y»AJ\L\ - oo) = (-ir/|Lp + O (l/|Lp) . (3.16) 

This estimate can be extended to any real value u > 1 of the index to. To make the content 
of Eq. (3.13) more transparent, we display in Fig. 1(a) the graphs of the analytic coupling 
for indices from —3 to and values of L in the range —3 to 3. Appealing to Eqs. (3.17) and 




FIG. 1: (a): Comparison of different curves for Ai, (L) as functions of the index v < 0, corresponding to 
various values of L, ranging from L = —3 to L = +3. The blue broken lines correspond to L < 0, whereas 
the red thick broken lines correspond to L > 0. The two dotted lines denote the results for L — ±1, with 
the line associated with the value L = 1 being closer to the zero line, (b): The same comparison for v > 
using the same line designations, but for values of L, ranging from L = —10 to L — +10 (dashed lines). The 
dotted lines here correspond to L — ±2, whereas the short-dashed lines represent the result for L = ±5. The 
blue solid lines in both panels show Av (0). 

(3.18) for negative values of the index u, one makes sure that, for L = 1 (red thick dotted 
line above the zero line), this function is equal to unity for all integer values of the index 
v = —m. On the other hand, for L < —1, the value of A_ m (L) depends on whether or 
not to is even or odd. For even values, it is positive, whereas for odd values it is negative, 
therefore giving rise to oscillations shown in Fig. 1(a). Note also that for values of L > 1, 
the oscillatory behavior of the graphs for A U (L) starts to be much less pronounced (red thick 
broken lines) because L m is positive for all positive values of to. The opposite behavior is 
exhibited for L < 0, as one sees from the blue broken lines. 

From Fig. 1(b), we observe that, in the region where A V {L) is smaller than unity (as 
explicitly indicated in the figure), this function is monotonic in v for v < 2. On the other 
hand, in the region where A V {L) > 1 — possible only for v < 1 and L < — this function 
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starts to be non-monotonic in u, so that there are two different points V\ and ^2, both 
corresponding to the same value of A V {V). Focusing on the values of A V {V) for L > 0, we 
see that all curves are monotonic in L and are bounded by an envelope represented by -4.^(0) 
(the blue thick solid line in Fig. 1(b)). If we consider only the interval of v G (0, 1), then the 
monotonicity property extends also to the negative values of L. 

Contrary to that case, the coupling A m (L) oscillates in L [4] for higher values of m > 2. 
These oscillations are not visible in Fig. 1(b) because of the smallness of the corresponding 
amplitudes. They appear due to rather general reasons: 

(i) the asymptotic conditions given by Eq. (3.16): A m (— oo) = A m (oo) = for m > 2; 

(ii) the differential relation between A m and A 1 , expressed in Eq. (3.1). 

Therefore, A m+2 has m zeros in the vicinity of the former "Landau pole" (L = 0) [4] — see 
Fig. 2. This property is rather unexpected from the point of view of standard power-series 
perturbation theory and will be discussed below in connection with Eq. (3.22). 
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FIG. 2: Comparison of different curves for A m (L) as functions of L, corresponding to various values of 
the index m, ranging from m — 2 to m — 5. To show all the details in a single plot, higher couplings are 
multiplied by numerical factors in order to normalize them to the scale of A^\l). The solid line shows 
A^\L), the blue dotted line corresponds to 8 • A^\l), the green short-dashed line to 60 • A^ (L), and the 
red long-dashed line to 480 • A^{L). 

This oscillation property of the coupling extends to A U (L) for all real values of the index 
v > 2. 

To reveal the relevance of this representation for physical applications, let us now consider 
A for some particular values of the index v. For the case of a negative index, the A- v play 
the role of the "inverse powers" of A\ that may be considered as the images of a~ u . Then, 
expression (3.13) can be rewritten in the form 

A-AL) = V - , (3.17) 

A- m (L) =hmA^ m+£ (L) = L m , at m = 0,l,2,... , (3.18) 

e— »0 

where we have taken into account that for v > 

F(z,u) = Li v (z), (3.19) 

with Li„ being the well-known polylogarithm function. It is worth remarking here that the 
"inverse powers" A- m (L) = L m coincide with the inverse powers of the original running 
coupling a~ m (L) = L m . 
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To make explicit the properties of Eq. (3.13), we convert this equation into a series 
representation, using Eq. (A3), to obtain 



F(z,l-") = I» 



In 



+ £)((!-„- r) !=M (3.20, 



for |ln(z)| < 2tt, where ((u) is the Riemann (^-function. Now we are in the position to 
express A v in the form of a series, i.e., 

1 °° (— LY 

A V {L) = -— C(l - v ~ rV—^- for \L\ < 2vr (3.21) 

because the "standard logarithms", contained in both parts of expression (3.13), mutually 
cancel, as one verifies by substituting Eq. (3.20) into Eq. (3.13). Then, we can state the 
following important corollaries: 

1. A u (0) = ~ ~p( ) " * s an en ^ re function of v. 

Of particular importance are the following values: 

AM = \ A 2 (0) = ^ .4,(0) =0, A(0) = -^, -4 5 (0) = (3.22) 

that coincide with the results provided in [4]. Note that *4.2n+i(0) = for n > 1 is due 
to the property ((— 2n) = 0, while the set of *4.2n(0) is alternating in sign [47]. These 
properties illustrate the details of the coupling oscillations in the vicinity of L = for 
index values m > 2. A convenient series representation of A m for an integer index m 
is presented and discussed in Appendix A (item 3). 

/■(I _ e _ r ) 

2. Taking into account the relation — lim — = 5 0r , (see, e.g., Eq. (A7)), one 

6-0 T(e) 

can take the limit 



A (L) = lim Au(L) = 1 , (3.23) 

v— »0 

dispensing with the constraint \L\ < 2n and proving assertion (2.8). 

3. Equation (3.18) can be re-derived from representation (3.21) in a way similar to that 
described in the previous item. 

The upshot of these considerations is that the linear space {A n } is now completed via the 
inclusion of the elements A v for any real values of the index v, so that one can take derivatives 
with respect to this continuous variable — dubbed 'index derivative'. 



C. Analytic images of products of coupling powers and logarithms 

To this point, we have considered only powers of the running coupling, adopting the 
viewpoint of the Shirkov-Solovtsov APT. Now we are going to consider more complicated 
expressions, like (a(i)) L m , where the power v is a real number and the power m is an 
integer, following the broader "analytization" principle of KS [38, 39]. 
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To compute this image, we have first to determine the image of a v \n(a), which can be 
rewritten as the derivative of a v with respect to v\ viz., 



a? ln(a) = — a u . 
du 



(3.24) 



Due to the linearity of the differential operator, this derivative can be directly applied to 
any element of the completed space {A n } to generate the corresponding image, \a v ln(a)] an , 
and define 



d 

dv 



def d 

du 



(3.25) 



In the following, we shall employ for the sake of simplicity a special notation for the deriva- 
tives with respect to the index of the non-power expansion and define 



d k 



dv k 



From Eqs. (3.13) and (3.25), we obtain 

Kln(«)L = - 



InL 



d 

dv 



F(e 



v] 



T(u) 



(3.26) 



(3.27) 



and taking multiple derivatives on both sides of Eq. (3.25), we compute the image of a u L m , 
like in Eq. (3.27). This procedure applies to any desired degree m of such terms. 

The extension to higher loops makes use of the APT expansion of higher-loop quantities 
in terms of one- loop ones. Before doing that, we consider first the image of ap) on the basis 
of the perturbation expansion, given in Eq. (B5) in conjunction with Eq. (3.26), to obtain 
the element A\ at the two-loop level: 



A? (L) = A? + Cl V A^l 2 + ci{V 2 + V 1 -l) A% + O (v 3 A?l 4 ) 
This formula can be readily generalized to any index v. 



(3.28) 



Af{L) = AW + Cl vVAV 1 + clv 
+ 0(V*A%) , 



V 2 + V 1 - 1 



A 



(i) 

u+2 



(3.29) 



where C\ = &i/6q is an auxiliary expansion parameter. The quality of the two-loop approxi- 
mation for the lowest index (cf. Eq. (3.28)) and higher indices (cf. Eq. (3.29)) will be analyzed 
numerically in the next section. Here it suffices to mention that the achieved accuracy is of 
the order of about 1% down to L — 0. 

To construct the image of (a(2)) L, cf. Eq. (2.10), we first perform the "analytization" 
of Eq. (B7) and then use Eq. (3.25) to arrive at the final expression 

l(a<n(L)y L] m = tftl(L) = 41 + cil>4 2) + O (a%) , (3.30) 

which can be recast, by means of the one-loop analytic coupling and with the aid of Eq. 
(3.27), in the form 



C^L)=A^ l -c 1 



(i) 



v 



h W-*M +Wl/)4 ., + ?n«-M-^) 



I» 



+ 0(t? 2 4 1 ] 1 ).(3.31) 
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The "analytization" of (a^)(L)Y ■ L 2 , expressed in terms of Eq. (B9), can be performed in 
an analogous way with the result 

C%{L) = A^ 2 + 2 Cl X?4 2 i + c\V 2 A^ - 2c\Af + O {vA%) . (3.32) 
The one-loop approximation of this two-loop expression is given by 

= 4% + Cl vVSSl + c\^±±V 2 A^ + O (VA%) . (3.33) 

A compilation of the required formulae to achieve the "analytization" of powers of the 
coupling in conjunction with logarithms at the two-loop order, is provided in Appendix B. 

Up to now we have studied expressions appearing in fixed-order perturbation theory. But 
similar considerations apply also to resummed perturbation theory. Indeed, first attempts 
to apply the "analytization" procedure of Shirkov-Solovtsov were already presented in [36, 
37, 39]. The crucial point here is how to deal with the requirement of analyticity when 
performing a Sudakov resummation. Because of the non-power series character of APT, 
resummation of (soft-gluon) logarithms does not lead to exponentiation. The latter can be 
retained only in the case of the so-called naive "analytization" [2], proposed in [36, 37]. The 
exact expression for the Sudakov factor is too complicated and too specific to be discussed 
in the present analysis. We, therefore, consider in Appendix C a simplified version of a 'toy 
Sudakov' factor that, nevertheless, bears the key characteristics pertaining to resummation 
under the assertion of analyticity. For clarity, we compare the basic ingredients of FAPT in 



TABLE I: Comparison of the standard PT, APT, and FAPT with L = In (Q 2 /A 2 ) . 



Theory 


Space 


Series expansion 


Inverse powers 


Multiplication 


Index derivative 


PT 




F(L) = £/ m e^(L) 

m 




/l v fl+U 

a (l) a (l) - a (l) 


_m i„fc 

a (0 m o (J ) 


APT 


urn 

l ) meN 


F( L ) = Y,fm Am (L) 

m 


No 


No 


No 


FAPT 


l J i/GR 


HL)=Y,fmA m ) {L) 

m 


A m m (L) = L m 


No 


V k A^ 



Table I with their counterparts in conventional perturbation theory and APT. More detailed 
expressions are shown in Table II in the next section. 



IV. VALIDATION OF THE NEW SCHEME 

A. Analytic verification of the one-loop spectral density 

An alternative way to derive Eq. (2.15) for the spectral density p u , is to compare two 
different representations for A u : one given by the dispersion relation, Eq. (2.3), and the 
other provided by the Laplace representation, Eq. (3.3). Then, we get 

A V {L) = H-^Tn do = r e ~ U Mt) dt . (4.1) 
Jo ° + v Jo 
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Next, we make a double Borel transformation of both representations, the Laplace one and 
that of the dispersion integral, the aim being to extract p v {(?)- This is done by applying first 
M 2 £?( M 2^q2) on both sides of Eq. (4.1) and then employing 



M 2 B 



(M 2 - 



<y + Q 2 



exp(- ( x/M 2 ), M 2 B (M 2^ Q2) 



A 2 ^ 

Q 2 



r(t) \m 2 ) ■ [ ^ Z) 



In the second step, we carry out one more Borel transformation, B(i/o--+i/m 2 ), to obtain 

'A 2 ysin(7rt) - 



A v {t) dt. 



(4.3) 



The final step is to substitute in Eq. (4.3) the expression for A v (t), given by Eq. (3.7), to 
arrive at the final result 



PM = 



1 



(n 2 + L 2 (a)) v/2 
1 



(tt 2 + LV)) 



v/2 



sin 



sin 



cr 



v arccos 



v arctan 



for L(a) > , 



7T 



L(a) 



(4.4) 
(4.5) 



where L(a) = In (a /A 2 ). To gain a more complete understanding of the role of the Landau 
pole remover in A v , it is important to remark that it does not contribute to the spectral 
density, the reason being that this part is not altering the nature of the discontinuity. The 
latter is solely determined by the term 1/L. One appreciates that expressions (4.5) and (2.15) 
coincide, as they should, hence establishing the equivalence between the two alternative 
extensions of the "analytization" procedure to fractional indices. The two-loop approximate 
expression for the spectral density is given in Appendix B. 



B. Verification of the two-loop approximations 

Now look specifically at the quality of the two-loop expansion in FAPT. In doing so, 
we define the following quantities with the help of an auxiliary parameter c\ and the index 
derivative V, (as in Eq. (3.28)): 



NLO, i.e., retaining terms of order c\ 

' ) • n DA 

NNLO, i.e., retaining terms up to order c\ 



AfPT(L) = 1 _^m±|M^) (4 . 6) 



A^\L) + Cl V A%( L) + cj (D 2 + P 1 - 1) A%{L) 

Af\L) 



For the corresponding quantities within the standard QCD perturbation theory, we use 
Eq. (B5) to obtain 
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. 1 


-Af APT (L) 


(a) ' 


0.08 






0.06 






0.04 






0.02 


\Af APT (L) 











L 



0.1 
.08 
.06 
. 04 
. 02 




(b) 



\ 



Af APT (L) 



01234567 



FIG. 3: (a): The dotted line corresponds to A FAPT (L) and the solid line to A FAPT (L). (b): The dashed 
line corresponds to A PT (L) and the solid line to A FAPT (L). 

• NLO, i.e., retaining terms of order c\ 



A-(L) 



^ a ( i) (L) + ci a 2 (1) (L) In a ( i) (L) 



a (2) (L) 



(4.8) 



NNLO, i.e., retaining terms up to order cf 



. PT 



O(i) (L) + ci o^j (L) In o(i) (L) + cf a 3 (1) (L) (in 2 a ( i) (L) + In a (1) (L) - l) 



(4.9) 

First, let us compare the transition from the NLO (cf. Eq. (4.6)) to the NNLO (cf. Eq. 
(4.7)) in FAPT (see Fig. 3(a)). One appreciates that by taking into account the NNLO terms, 
a significant improvement of the convergence quality of the FAPT series is achieved. Indeed, 
even at Q 2 = A 2 , which corresponds to L — 0, the error of truncating the FAPT series at 
the NLO is about 5%, while by taking into account the NNLO correction this error becomes 



0.06- 



0.04 



0.02 - 





(a) - 




\ X 
\ \ 
\ \^ 

\ >v 

\ X S( ^ 







2 L 4 



10 




f2) f2VFA.PT 

FIG. 4: (a) Comparison of different results for A\ \L). The solid line corresponds to A 2 (L), com- 



puted analytically via Eq. (3.29); A 



(2) ;num 



(L) (dashed line) is derived by means of a numerical integration. 



The dotted line represents the available results of the numerical procedure of Magradze in [23] . (b) Compar- 
ison of FAPT and standard QCD PT with respect to the fractional index (power) of the coupling. The solid 
line represents ^4Q 2 g2 FAPT (L), computed analytically via Eq. (3.29), while the dashed line stands for a° 2 ^ 2 (L). 
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even smaller than 0.5%. In Fig. 3(b) we show the relative quality of these approximations 
concerning the loop expansion between the standard perturbation theory and FAPT, as 
quantified by Eqs. (4.7) and (4.9). One appreciates the strong suppression of A3 APT (L) 
relative to its conventional analogue in the small L region, say, below approximately L = 2. 

The same comparison can be realized for A^\ using Eq. (3.29). Indeed, we demonstrate in 
Fig. 4(a) the quality of this FAPT expansion in comparison with the results of the numerical 
integration of the NLO spectral density p2 (for more details, we refer to Appendix B and 
[5]) in the dispersion-integral representation, provided by Eq. (2.3). In this graphics, we 
also display the results obtained numerically by Magradze in [23]. The message from Fig. 
4(a) is quite clear. Our analytic (solid line) and our numerical calculation (dashed line) 
are in mutual support, while the results of [23] differ considerably with respect to both the 
magnitude and the trend of the negative values of L. The good convergence of the proposed 



o. 



o. 



o . 



1 




8 












6 












4 












2 



LdSc5 








L 




12 


3 


4 


5 


6 7 



FIG. 5: Analytic images Al (L) at the two-loop level for some characteristic index values, calculated with 
NNLO FAPT. The solid line corresponds to v = 0, the dash-dotted one to v = 1, and the lowest line to 
v = 2. All other lines correspond to v = k/A with k = 1, 7. 

a We wish to thank D. V. Shirkov for suggesting to us this form of presentation. 

series for Af\ (Eq. (3.29)), that had been demonstrated above, can be traced to the basis of 

APT. Indeed, this non-power expansion of the quantities W[l in terms of Am has a finite 
radius of convergence, the reasons being discussed in Ref. [49]. 

To give the reader an impression of the dependence on L of Wi\L), we show in Fig. 4(b) 
a comparison of this quantity with its counterpart in standard QCD, namely, [a(2)(L)l • For 
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TABLE II: Calculational rules for FAPT with L = In (Q 2 /A 2 ) , m G N, and v G E. 



Standard QCD PT 


QCD FAPT 




-<(£)=!- l' , 
L e n — I 




1 J " L" r(i/) 


afo(L) ln m [a (0 (L)] 






A?\L) =AP(L) + c 1 vVAU 1 (L) + cIv \^±^V 2 + V- ij^L) 

+ ^ 3 4 1 ] 3 ) a 


a^ 2) (L)L 


4 2 J (l) = 4 2 i(D + d p4 2) (l) + o (v 2 4 2 i) 6 

« ^ w - - [ ln(L) - * w + *wv\l) + V XA " 1 


a^ 2) (L)L 2 


C^(L)=A?l 2 (L) + 2c 1 VA?l 1 (L)+c 2 1 V 2 A^(L) - 2c 2 A^(L) 
+ 0{VA {2) +I ) b 

kA^ + cvva^ + c 2 [ 2 ~ u 2 +A ] v 2 a£Xl) 


exp [— xa(L)] 


e -*/L + e ~ mL Jl (2 ^ } for L > c 



a Note that in evaluating this expression in the next line we use Eqs. (B6) and (3.13). 
fc Note that in evaluating this expression in the next line we use Eqs. (B8), (B9), and (3.13). 
Tor the derivation of this expression, we refer to Appendix C. 



the purpose of illustration, we select the value v — 72/ (2&o) ~ 0.62, which corresponds to the 
1-loop evolution exponent of the non-singlet quark operator of index 2, entering a number 
of applications in DIS and also various exclusive reactions [43, 50, 51]. 

In support of our two- loop approximation (within FAPT), we display in Fig. 5 results for 

(2) 

the analytic images Ah with v — k/A and k — 0, 1, . . . , 8. We observe the same monotonic 
pattern, i.e., no crossing, of curves, found already for the one-loop case, shown in Fig. 1(b). 
We investigated numerically the range of negative values of L and found that crossing appears 
only for indices v < 1, again in close analogy to the one- loop case. 

The pivotal results of this paper are collected in Table II, where we provide the reader with 
explicit calculational rules to connect the standard QCD perturbation theory with FAPT. 
We stress that the presented algorithm has broad applications in phenomenology and can 
play a major role in the perturbative analysis of observables both in inclusive and exclusive 
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QCD reactions. 



V. CONCLUSIONS 

With hindsight we can say that the requirement of analyticity at the amplitude level of 
hadronic quantities in QCD, expressed by Karanikas and Stefanis [38], is instrumental in 
improving perturbation-theory calculations. First, as shown in an accompanying paper by 
two of us (A.P.B. and N.G.S.) together with Karanikas [40], it enables one to minimize both 
the sensitivity on the renormalization scheme and scale setting and also the dependence 
on the factorization scale. The reason for this latter advantage is that it includes into the 
"analytization" procedure not only the powers of the running coupling, but also logarithms 
(or exponentials) that may contain the momentum scale with respect to which analyticity is 
required. Second, starting at this point, we have shown in the present work that invoking this 
analyticity principle gives rise to a generalization of the original APT to fractional powers 
of the coupling. As a bonus, this approach improves the convergence of the perturbative 
expansion significantly — see Fig. 3. 

Our main goal in this analysis was to work out in mathematical detail the procedure for 
determining analytic expressions for any real power of the running coupling and delineate 
the main results. To keep our presentation as general as possible, we have purposefully 
refrained from considering specific examples and concentrated instead on generic features 
and expressions. In this vein, we have discussed products of the running coupling (or its 
powers) with (powers of) logarithms, which are typical for contributions encountered in 
higher-order corrections of QCD perturbation theory, or when taking into account evolution 
effects via the renormalization-group equation. Similar logarithmic terms also appear in 
calculations employing light cone sum rules [42, 43]. In an analogous way, we have discussed 
the resummation of non-power series in the analytic images in order to capture the key 
features of Sudakov resummation of soft-gluon effects (see Appendix C). All these elements 
of FATP, required for further applications of this formalism to improve the calculation of 
any hadronic amplitude at the two-loop level, are collected in Table II. In this context, we 
mention that we have developed approximate expressions for the two-loop analytic images 
in terms of one-loop quantities that can facilitate practical computations significantly. 

In conclusion, this report has emphasized rigorous methods rather than specific applica- 
tions. A first example of the present framework is discussed in [40], focusing on the topic 
of the renormalization-scheme and factorization-scale independence of the electromagnetic 
pion form factor, relative to its treatment within the standard QCD perturbation theory or 
original APT. The methods presented here are intended to be used in the low-to-medium 
momentum range where the standard perturbative approach faces the problem of the Landau 
pole and in processes or under circumstances where the original APT is insufficient because 
it is tied to integer powers of the coupling. We believe that our assortment of analytic ex- 
pressions for a variety of expressions ranging from any real powers of the coupling to more 
complicated products containing logarithms, provides sufficient evidence for the usefulness 
of the approach for higher-order perturbative calculations. 

Acknowledgments 

We wish to thank D. V. Shirkov for valuable discussions and comments. Two of us (A.P.B. 
and S.V.M.) are indebted to Prof. Klaus Goeke for the warm hospitality at Bochum Univer- 



17 



sity, where the major part of this investigation was carried out. This work was supported in 
part by the Deutsche Forschungsgemeinschaft, the Heisenberg-Landau Programme (grant 
2005), and the Russian Foundation for Fundamental Research (grants No. 03-02-16816, 03- 
02-04022 and 05-01-00992). 



APPENDIX A: ANALYTIC PROPERTIES OF -u, 1) AND A v 

1. The function $(z, — is, i) can be determined by means of the analytic continuation of 
the series [47] 



$>(z,s,i) = 



m=0 



(m + i) s 



for \z\ < 1, s ^ 1,2,3, ... 



(Al) 



in both variables z and s. This analytic continuation for every fixed s, which is not a positive 
integer, determines <E> as an analytic function of z, regular in the plane with a cut along the 
axis (1, oo), and for every fixed z in the cut plane as an analytic function of s being regular, 
except possibly at the points s = 1, 2, 3, . . . , as was mentioned in [47]. We can improve 
this statement for s > 1 using Eq. (3.19) to obtain 



$(z,s,i) = - 



Hz,s,l)-J2 — 

m s 



m>l 



1 



i-1 



m>l 



Li s (z)-Y — 



(A2) 



One appreciates that there are no singularities for any positive integer values of s. Hence, 
we can conclude that $(z, s, 1) is an analytic function in s at any fixed z on the cut plane. 
Moreover, the function A U (L), see expression (3.11), 



ML) = JZ - 



$(e- L ,l-M) 

r(iz) 



has no poles in v and is, therefore, an entire function in v. 

2. There is a useful series representation for s, 1) [47] (cf. (3.20)); viz., 



$(z,s,l) 



1 

z 



r(i- s ) 



In 



-I s-l 



ln r (^) 



r=0 



T ! 



(A3) 



for |ln(z)| < 27r, s ^ 1,2,3,... that allows one to continue $(z,s,l) for integer positive 
s = m values by means of the limit s = m + e, e — > in Eq. (A3). To take this limit, 
one should expand in e the first term in the square brackets in expression (A3), which is 
proportional to T(l — m — e). The other singular term appears in the sum and is proportional 
to C(l + e) — 1/e — ip(l) + 0{e). The singularities, contained in both these parts, mutually 
cancel. 

3. The expansion of A V {L), Eq. (3.21), is simplified for an integer index v = m > 1 to 
read 



A m {L) 



' r=0 



B 



m+r 



(m — 1)! ^ (m + r)r\ 



-L) r for \L\ < 2tt, 



(A4) 



where B m are the Bernoulli numbers. From the property B 2n +i = it follows that A2n(L) 
is an even function of its argument, while A2 n +i(L) is an odd one. 
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Note here that the pole remover F(e~ L , 1 — m) in expression (3.13) reduces to elementary 
functions for the case of an integer index. Indeed, according to Eq. (3.1), the operator —d/dL 
shifts the second argument of the function F by unity, i.e., m — > m + 1, to get 

Taking into account that F(e~ L , 0) = (e L — l)^ 1 and applying the previous relation m times 
we arrive at 

^ i --™»-(-^)'" F(e " 0) K-^)'"^T- (A5) 

This representation leads to an exponentially suppressed asymptotic limit for the function 
F(e~ L , — m); viz., 

F(e- L ,-m) ~e~ |L| . (A6) 

4. We supply here the Lindelof formula [52] 

, 1 1 /"°° sin [z/arctan(t)l (it . . „, 

CM ^ 2 + — + X (i + - 1) ' (A7) 

which fixes ((u) as an analytic function with a simple pole at v — 1. This representation for 
C(z/) has been used in Sec. Ill B. 

5. Now we are in the position to supply also the analytic images of the coupling in the 
timelike regime for L(s) = log(s/A 2 ) > 0, employing the notation of [6, 53]: 



DC 



da 1 [°° sin [z/arctan (7r/L)] 

° * Jl(s) (tt 2 + L 2 ) 



M*)= I -pM = -\ dL - - — ^ » . (A8) 



This integral can be evaluated to provide a result analogous to A V (L) for the spacelike 
regime; namely, 

Ms) = ^M^-^arctanWL^))] 

7T(Z/-1)(7T 2 + L( S ) 2 ) ( ^ 1)/2 

A similar expression for the timelike coupling has been obtained before in [54] using the 
"contour-improved resummation technique" . 



APPENDIX B: "ANALYTIZATION" OF POWERS OF THE COUPLING MUL- 
TIPLIED BY LOGARITHMS 

1. The expansion of the /3-function in the NLO approximation is given by 

* j«jm 2 _j?M\\ (B1) 



dL \ An J \ An J \ An 

where L = ln(/z 2 /A 2 ) and 



11 4 34 / 20 \ 

b = — C A - - T R N f , h = — Ci - UC F + — C A \ T R N f (B2) 
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with C F = (iV c 2 - 1) /2N C = 4/3, C A = N c = 3, T R = 1/2, and N f denoting the number of 
flavors. Then, the corresponding two- loop equation for our coupling a = bo a /(Air) looks like 



da 



(2) 



bi 



dL 



= -af 2) (L) [l + ci a(2)(L)] with ci _ 2 . 



The renormalization-group solution of this equation assumes the form 



a (2) 



+ Ci In 



1(2) 



= L . 



(B3) 



(B4) 



order a 3 



1 + CiO( 2 )_ 

Then, for the expansion of ap) (L) in terms of (L) = 1/L we have, retaining terms of the 

3 

(1)' 

0(2) = O(i) + ci a 2 {l) lna ( i) + c 2 a 3 1} (in 2 a ( i) + lna (1) - l) + C(a( 1} ln 3 (ai)) . (B5) 
2. Now, for the product [a(2)]"l/, we obtain from (B4) 

°(2) 



{ a (2)Y L = (a(2))" 1 + (a(2))"ci In 



1 + cia (2 ) 



(B6) 



Expanding the logarithmic term ln[l+cia( 2 )], while retaining terms of order aj^ , ln(a(2)); 
viz., 



(a(2))"L = (0(2))" + C! (a {2 )Y H a (2)) - Oia^ 1 ) 
and, finally, expanding the coupling a( 2 ) in terms a = am, we find 



(a m ) v L 



+ u a u Cl ln(a) + 0{a v+1 \n 2 (a)) 



(B7) 



(B8) 



Calculating (a^)Y L 2 in an analogous way, we derive 

(a (2 )Y L 2 = i a (2)Y~ 2 [l + cia (2) In (a (2) )] 2 - 2 c 2 a^ 2) - C(a^+ 1 ln(a (2) )) 



= a 



.f-2 



+ z/ a 1 ' ci ln(a) + 



z/ + 4 



a u c{ ln 2 (a) + 0(a" ln(a)) . (B9) 



3. We consider here the spectral density p v {p) beyond the leading-order approximation. 
At the /-loop level, pv\o) can always be presented in the same form as for the leading order, 
given in Eq. (2.15), i.e. 



— Im \c 

7T L 



(BIO) 



but keeping in mind that the phase <f^ and the radial part have a multi-loop content. 
At the two-loop level, one should, strictly speaking, deal with the imaginary part of the 
Lambert function W-\ (see [20]) because the exact solution of Eq. (B4) can be realized in 
terms of the Lambert function. Instead of following this procedure, we can alternatively take 
the well-known first-iteration solution of Eq. (B4) that provides us with sufficient accuracy 
the following result: 



a {2) (L) a^(L) 



L + ci ]n(L + d) 



(Bll) 
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For this approximate solution a^f, we have 

{R { 2)(v)) 2 = (L{a) + d In (y<M*) + Cl ) 2 + tt 2 )) 2 + (tt + Cl 0(L(<r))) 2 , (B12) 

L(cr) + ci In (v/(L(a) + Cl ) 2 + tt 2 



</?( 2 ) ^crj = arccos 
<f)(L(cr)) = arccos 



L(a) + ci 



(B13) 
(B14) 



y(L(a) + ci) 2 + 7r 2 

with L(a) = In (a/A 2 ). The spectral density p^=^ a (a) with the phase and the radial part 

(2) 

from Eqs. (B12)-(B14) appears to be very close to the numerical, but exact result for p\ (a), 
based on W-\ — see, e.g., [5]. 



APPENDIX C: "ANALYTIZATION" OF THE TOY MODEL FOR SUDAKOV 
RESUMMATION 

Here we discuss the analytic image of expression F$(x, L) = exp [-ia(i)(I)] , which orig- 
inates as a part of the procedure of the Sudakov resummation, where x is a free parameter. 
We consider the following example, serving as a "toy model" for this resummation: 



{F S (x,L)}, 



[-xa {l) (L)]'' 



m=l 



ml 



-X] 



m=l 



ml 



(Cl) 



One can verify that for the asymptotic limits of L, Eq. (Cl) reduces to the evident forms: 



l + X)(-aO 



m=l 



ml 



1 — x for L — > — oo ; 
1 for L-> +oo . 



(C2) 
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Ft\x,L) 


(a) 






x = L0. _ - - ' " ' 

~ — 




y 


x = — ^ / 
~~ ~~ / 
x = 2.0- ^ 


L - 



-10 



. 8 
. 6 
. 4 
0.2 

10 -io 



F s (x = l,L) 


(b) 




/ L 



10 



FIG. 6: (a) Functions Fg an ^ (x, L) calculated for different values of x using the one-loop FAPT. The horizontal 
solid line corresponds to x = 0, while the broken lines represent the specific x values shown above them, 
(b) Comparison of the toy Sudakov function, calculated within the present framework (solid line) , with the 
result of the naive "analytization" [2, 36, 37] (dashed line), and that from conventional QCD perturbation 
theory (dotted line). Note the different scale for the ordinate relative to panel (a). 
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The first asymptote on the RHS of Eq. (C2) appears due to the equality Am(— oo) = Si m , 
see, for instance, [4]. The second asymptote is due to the property that in the UV regime 
APT reduces to the standard perturbation theory. Both properties are illustrated in Fig. 
6(a). In Fig. 6(b) we compare different versions of the toy Sudakov function, obtained within 
the present framework (solid line), using the naive "analytization" [2, 36, 37] (dashed line), 
and conventional QCD perturbation theory (dotted line). 

On the other hand, restricting the loop order to I — 1, one can derive an explicit expression 
for Eq. (CI) in the region of L > 0_, which is based on the Laplace representation given by 
Eq. (3.3) and expression (3.8) for Ai, namely, 

Ft\x, L) = e~ x ' L + ^J2 e ' Lm ^^f^ . L>0 - (C3) 

171=1 V 

The perturbative part of A m in (CI) reproduces exactly the asymptotic expression 
exp [— x a(i)(L)] on the RHS of Eq. (C3), while the pole remover generates the sum of the 
exponents in —L, weighted by the Bessel functions, Ji, exhibiting how the large L behavior 
is violated. 



S. J. Brodsky, G. P. Lepage, and P. B. Mackenzie, Phys. Rev. D28, 228 (1983). 

A. P. Bakulev, K. Passek-Kumericki, W. Schroers, and N. G. Stefanis, Phys. Rev. D70, 033014 

(2004); [Erratum-ibid. D70, 079906 (2004)]. 

D. V. Shirkov and I. L. Solovtsov, Phys. Rev. Lett. 79, 1209 (1997). 
D. V. Shirkov, Theor. Math. Phys. 119, 438 (1999). 
I. L. Solovtsov and D. V. Shirkov, Theor. Math. Phys. 120, 1220 (1999). 
D. V. Shirkov, Eur. Phys. J. C22, 331 (2001). 

D. V. Shirkov and I. L. Solovtsov, Phys. Part. Nucl. 32S1, 48 (2001). 
G. Parisi and R. Petronzio, Nucl. Phys. B154, 427 (1979). 
J. M. Cornwall, Phys. Rev. D26, 1453 (1982). 

M. B. Gay Ducati, F. Halzen, and A. A. Natale, Phys. Rev. D48, 2324 (1993). 
A. C. Mattingly and P. M. Stevenson, Phys. Rev. Lett. 69, 1320 (1992). 
A. C. Mattingly and P. M. Stevenson, Phys. Rev. D49, 437 (1994). 
D. V. Shirkov, hep-ph/0003242 (unpublished). 
D. V. Shirkov, hep-ph/0009106 (unpublished). 

D. V. Shirkov, hep-ph/0408272 (unpublished). 

I. L. Solovtsov and D. V. Shirkov, Phys. Lett. B442, 344 (1998). 

A. P. Bakulev, A. V. Radyushkin, and N. G. Stefanis, Phys. Rev. D62, 113001 (2000). 

G. Grunberg, hep-ph/9705290 (unpublished). 

E. Gardi, G. Grunberg, and M. Karliner, JHEP 07, 007 (1998). 

B. A. Magradze, Int. J. Mod. Phys. A15, 2715 (2000). 

D. S. Kourashev and B. A. Magradze, hep-ph/0104142 (unpublished). 

D. S. Kurashev and B. A. Magradze, Theor. Math. Phys. 135, 531 (2003). 

B. A. Magradze, hep-ph/0305020 (unpublished). 

A. V. Radyushkin, JINR Rapid Commun. 78, 96 (1996). 

N. V. Krasnikov and A. A. Pivovarov, Phys. Lett. B116, 168 (1982). 

H. Contopanagos and G. Sterman, Nucl. Phys. B419, 77 (1994). 
M. Beneke and V. M. Braun, Phys. Lett. B348, 513 (1995). 



22 



[28] P. Ball, M. Beneke, and V. M. Braun, Nucl. Phys. B452, 563 (1995). 

[29] Y. L. Dokshitzer, G. Marchesini, and B. R. Webber, Nucl. Phys. B469, 93 (1996). 

[30] S. Catani, M. L. Mangano, P. Nason, and L. Trentadue, Nucl. Phys. B478, 273 (1996). 

[31] L. Magnea, Nucl. Phys. B593, 269 (2001). 

[32] E. Gardi, Nucl. Phys. B622, 365 (2002). 

[33] A. V. Nesterenko, Int. J. Mod. Phys. A18, 5475 (2003). 

[34] A. V. Nesterenko and J. Papavassiliou, Phys. Rev. D71, 016009 (2005). 

[35] A. C. Aguilar, A. V. Nesterenko, and J. Papavassiliou, J. Phys. G31, 997 (2005). 

[36] N. G. Stefanis, W. Schroers, and H.-C. Kim, Phys. Lett. B449, 299 (1999). 

[37] N. G. Stefanis, W. Schroers, and H.-C. Kim, Eur. Phys. J. C18, 137 (2000). 

[38] A. I. Karanikas and N. G. Stefanis, Phys. Lett. B504, 225 (2001). 

[39] N. G. Stefanis, Lect. Notes Phys. 616, 153 (2003). 

[40] A. P. Bakulev, A. I. Karanikas, and N. G. Stefanis, Phys. Rev. D72, 074015 (2005). 

[41] K. A. Milton and I. L. Solovtsov, Phys. Rev. D55, 5295 (1997). 

[42] A. Schmedding and O. Yakovlev, Phys. Rev. D62, 116002 (2000). 

[43] A. P. Bakulev, S. V. Mikhailov, and N. G. Stefanis, Phys. Rev. D67, 074012 (2003). 

[44] N. N. Bogolyubov, A. A. Logunov, and D. V. Shirkov, Soviet Physics JETP 37, 574 (1960). 

[45] R. Oehme, Phys. Lett. B252, 641 (1990). 

[46] B. A. Magradze, hep-ph/0010070 (unpublished). 

[47] H. Batemann and A. Erdelyi, Higher transcendental functions (Batemann Manuscript Project) 

(McGraw-Hill, New York, 1953). 
[48] F. Olver, Asymptotics and special functions (Academic Press, New York, 1974). 
[49] S. V. Mikhailov, hep-ph/0411397 (unpublished). 

[50] A. P. Bakulev, S. V. Mikhailov, and N. G. Stefanis, Phys. Lett. 578, 91 (2004). 
[51] A. P. Bakulev, S. V. Mikhailov, and N. G. Stefanis, Phys. Lett. B508, 279 (2001). 
[52] E. T. Whittaker and G. N. Watson, A course of modern analysis (Uni. Press, Cambridge, 
1927). 

[53] D. V. Shirkov, Theor. Math. Phys. 127, 409 (2001). 

[54] D. J. Broadhurst, A. L. Kataev, and C. J. Maxwell, Nucl. Phys. B592, 247 (2001). 



23 



